Enhanced Stability of Superheavy Nuclei due to High-Spin Isomerism 
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Configuration-constrained calculations of potential-energy surfaces in even-even superheavy nu- 
clei reveal systematically the existence at low excitation energies of multi-quasiparticle states with 
deformed axially symmetric shapes and large angular momenta. These results indicate the preva- 
lence of long-lived, multi-quasiparticle isomers. In a quantal system, the ground state is usually 
more stable than the excited states. In contrast, in superheavy nuclei the multi-qausiparticle exci- 
tations decrease the probability for both fission and a decay, implying enhanced stability. Hence, 
the systematic occurrence of multi-qausiparticle isomers may become crucial for future production 
and study of even heavier nuclei. The energies of multi-quasiparticle states and their a decays are 
calculated and compared to available data. 
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It has long been a fundamental question as to what are 
the maximum charge and mass that a nucleus may at- 
tain. Recently, significant progress has been made exper- 
imentally in the synthesis of the heaviest elements 0, • 
According to classical physics, elements with Z > 104 
should not exist due to the large Coulomb repulsive force. 
The occurrence of superheavy elements with Z > 104 is 
entirely due to quantal shell effects. Theoretically, pre- 
dictions have been made that there should be an island 
of stability of superheavy nuclei around N — 184 and 
Z = 114, 120 or 126, depending on the model employed 



Information concerning the excited-state structure of 
superheavy nuclei is scarce, so that the testing of struc- 
ture calculations has been limited. A promising area of 
progress is represented by the observation of rotational 
7-ray transitions in 2 q 2 ' 2 4 No H, @ and 2 gg' 256 Fm [To| . 
showing the deformed character of these nuclei, which 
can be described theoretically [ijj . The 7-ray transitions 
from non-collective one-quasiparticle states form another 
important probe into the structure of superheavy nuclei 
[I2j. By measuring the properties and decays of quasi- 
particle states, essential information about the nucleonic 
orbits close to the Fermi surfaces can be gained. Indeed, 
detailed knowledge of single-particle states is required for 
the development of reliable nuclear structure models. 

In the present work we consider multi-quasiparticle 
(multi-qp) states formed by breaking pairs of nuclc- 
ons. Unpaired nucleons couple their angular momenta 
to states with total spin projection A onto the symme- 
try axis in deformed, axially symmetric nuclei. The A 
values can be large, and, due to the high degree of their 
A forbiddenness, 7-ray transitions from high- if into low- 
A states are strongly hindered, leading to the formation 
of long-lived high-A isomers. In addition, there are in- 
fluences on fission and a-particle emission, which can be 
important in increasing the survival probability of su- 



perheavy nuclei. Indeed, high-A isomers might become 
crucial in the exploration of even heavier nuclei. In the 
present work, we focus on high-spin isomerism and the 
angular-momentum effect on the stability of even-even 
superheavy nuclei. 

Configuration-constrained potential-energy-surface 
(PES) calculations [l^J have been performed to de- 
termine the deformations and excitation energies of 
multi-qp states. Single-particle levels are obtained from 
the non-axial deformed Woods-Saxon potential with 
the set of universal parameters [14J. In order to reduce 
the unphysical fluctuation of the weakened pairing field 
(due to the blocking effect of unpaired nucleons) an 
approximate particle-number projection has been used 
by means of the Lipkin-Nogami method |Tfij with pai ring 
strengths determined by the average gap method [lfij . 
In the configuration-constrained PES calculation, it is 
required to adiabatically block the unpaired nucleon 
orbits which specify a given configuration. This has 
been achieved by calculating and identifying the average 
Nilsson quantum numbers for every orbit involved in a 
configuration The total energy of a state consists 

of a macroscopic part that is obtained with the standard 
liquid-drop model ^3 an d a microscopic part which is 
calculated by the Strutinsky shell-correction approach, 
including blocking effects. The PES is calculated in 
the space of quadrupole (/?2, 7) and hexadecapole 
(Pi) deformations. The configuration-constrained PES 
calculation can properly treat the shape polarization 
due to unpaired nucleons. 

Following predictions from theory, many experiments 
[8J, 13, LLfJ have demonstrated the systematic existence of 
deformed superheavy nuclei. In fact, almost all the su- 
perheavy nuclei found experimentally are believed to be 
deformed, and the stability of superheavy nuclei is en- 
hanced due to deformed shell effects. Furthermore, de- 
formed high- A orbits close to the Fermi surfaces can form 
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FIG. 1: Calculated excitation energies for two-quasineutron 
states. 
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FIG. 2: Calculated excitation energies for two-quasiproton 
states. 



various high-i-T multi-qp states at low excitation ener- 
gies. In the present work, we are interested in possible 
isomeric states. The combination of high K, low energy 
and axially deformed shape provides the necessary condi- 
tions for the formation of isomers. Additionally, a parity 
opposite to that of the ground state (g.s.) can further 
reduce the transition rate from an excited state to the 
g.s., or to a member of the g.s. band (g.s.b.). Figs. 
and display our calculations of low-energy high-X two- 
qp states with negative parity for even-even Z > 100 
nuclei. These investigated states have deformed, axially 
symmetric equilibrium shapes. 

In a recent experiment [Tsj . a decays from isomeric 
states in fipDs were observed with a half-life of 6.0IJ' 2 , 
ms that is significantly longer than the IOO^q fis of the 
corresponding g.s. Possible configurations for these iso- 
mers were proposed to be K r = 9" (z/^-"[725]®| + [613]) 

and/or 10~ (z/^~[725] ® | + [615]) [H|. Our calculations 
confirm these assignments and indeed show that the 9 _ 




FIG. 3: Schematic figure for the a decays corresponding to 
the reported decay events [Tsf from 270 Ds to 266 Hs. Differ- 
ent options for the initial isomer are indicated. The energies 
(in MeV) are not drawn to scale. The configurations are: 
u9- = z/^-~[725] g> | + [613]; vl0~ = u^~[725] ® | + [615]; 
7rl0~ = 7r^- + [615] <g> I [5051- The values in parentheses are 
the experimental energies |18| . 



and 10~ two-quasineutron states exist systematically at 
low excitation energies (see Fig. ^| with axially symmet- 
ric deformations around (3% ~ 0.2. In 270 Ds, the calcu- 
lated energies of the 9~ and 10~ isomers are 1.39 and 

I. 26 MeV, respectivel y, ag reeing with the experimental 
estimate of 1.13 MeV [lg- Furthermore, configuration- 
constrained PES calculations allow us to determine the 
Q a value of the a decay from a given configuration of the 
parent nucleus to a given state of the daughter. Then, the 
a-particle kinetic energy is obtained from E a = Ap A 4 Q a 

allowing for the recoil correction (where A p is the mass 
number of the parent nucleus). Fig. [3] shows schemat- 
ically the a decays from the g.s. and isomers in 270 Ds. 
The calculated a decays of the two-quasineutron z/9~ and 
i^lO isomers to the corresponding isomers of 266 Hs have 
very similar a-particle energies, giving an average value 
of 11.06 MeV. This is in good agreement with the ob- 
served energy of 10.95 MeV |l|. The u9~ and ul0~ 
isomeric a decays to the g.s. of 266 Hs also have similar 
a energies with an average value of 11.93 MeV, in agree- 
ment with the measured energy of 12.15 MeV [l8]. One 
also expects that the isomer can feed by a-emission into 
other states with low angular momenta. 

Also from 270 Ds, an observed a-emission with E a = 

II. 15 MeV was found in coincidence with a 218 keV 7- 
ray signal, and was suggested to decay to the 8 + member 
of the rotational g.s.b. of 266 Hs 18] (the 7-ray energy 
was found to be close to the calculated energy of the 
transition from the 8 + to 6 + states of the g.s.b. in 266 Hs). 
We calculate that the 8 + to 6 + transition has an energy 



3 



of 234 keV. The calculated a energies of the v9 and 
i^lCT decays to the 8+ state are 11.42 and 11.29 MeV, 
respectively. In 270 Ds, however, we found another low- 
lying 10 _ state that has a two-quasiproton configuration 
7r "T + [^15] ® | [505] with a calculated excitation energy 
of 1.07 MeV, giving E a {nl0- -> 8+) = 11.10 MeV, in 
good agreement with the observed value of 11.15 MeV. 
In 266 Hs, the 7rl0~ state has a relatively high excitation 
energy of 2.54 MeV. 

To understand the above observations, we note that 
the probability for a decay is determined by two factors: 
i) an exponential dependence on the height of the bar- 
rier that the a particle needs to penetrate and ii) the 
preformation factor, characterising the ease with which 
an a particle is formed at the surface. The latter is pro- 
portional to the pair density at the Fermi surface |lflj. 
implying a considerable reduction for two-qp states, as 
is the case for high-if isomers. An a decay between 
two different isomeric configurations is in principle for- 
bidden. Furthermore, if the a particle carries a certain 
angular momentum, the barrier is increased due to the 
centrifugal potential. The increase of the barrier height 
and width reduces the a-emission probability and hence 
increases the lifetime. However, due to the exceptionally 
high energy of the a particle for the case of superheavy 
nuclei, the centrifugal barrier is of less importance, and 
the decay from an isomer to the g.s. can in fact com- 
pete with the decay to the corresponding isomer. This is 
consistent with the observation of the relatively long life- 
times for the a decays of the isomers in 270 Ds and 266 Hs 
|l8j. Therefore, the high- if isomerism can enhance the 
stability of such nuclei against a-emission. 

High- if isomers have also been found in 250 ' 256 Fm and 
254 No. In 250 Fm, an 1.8-s isomer [2ij was observed at an 
excitation energ y o f 1.0 MeV [13 . The 254 No nucleus has 
a 0.28-s isomer :2Q] with an energy > 0.5 MeV [nj, and 
there is recent evidence for associated collective struc- 
ture plj . However, experiments have not determined 
spins and parities for these two isomers. A likely con- 
figuration in 250 Fm is 7r| + [633] (g> ? [514] resulting in a 
7~ state. As shown in Fig. |2 the calculated excitation 
energy of the 7~ configuration is 1.01 MeV, which agrees 
well with the experimental value of 1.0 MeV ^3- A 
two-quasineutron 8~ configuration, ^| + [624] ® | [734], 
is also possible with a calculated energy of 0.97 MeV 
in 250 Fm. Our calculations show that 8~ states exist 
systematically in the N = 150 isotones with excitation 
energies around 1.0 MeV. In 254 No, alternative configura- 
tions for the observed isomer are ^| + [613] ® | [734] and 
tt|"[514] ® f + [624] with = 1.12 and 1.48 MeV, re- 
spectively. The low- lying i/| + [613] <& § [734] states exist 
systematically in the N — 152 isotones with calculated 
energies around 1.1 MeV. 

Another interesting example is the 256 Fm nucleus in 
which a 70-ns, 7~ isomer was observed to undergo spon- 
taneous fission [22]. The 7~ isomer was populated from 
the j3 decay of an 8 + isomer of odd-odd 256 Es. The 
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FIG. 4: Calculated potential energies versus deformation 
in 256 Fm. The solid and dashed lines represent the ground 
state and the 7r| + [633] ® | [514] isomer, respectively. At 
each fa point, the energy has been minimized with respect to 
the 7 and /?4 deformations. 

256 Es 8 + isomer has a half-life of 7.6 hours, which is 
much longer than the 25 minutes of the corresponding 
1 + g.s. We calculate that the excitation energy of the 7~ 
(tt| + [633] <g> |~[514]) isomer in 256 Fm is 1.1 MeV (see 
Fig. |2J) about 0.3 MeV less than the experimental value 
of 1.42 MeV 0. 

In 256 Fm, the observed partial fission half-life of the 
7~ isomer is 0.8 ms, which is remarkably longer than 
the expected value of 2.5 /is assuming no hindrance from 
unpaired nucleons |22|. In Ref. j2^|, it has been shown 
that the spontaneous-fission half-life increases consider- 
ably due to the effect of the unpaired nucleon in odd-mass 
nuclei. This can be understood in terms of reduced su- 
perfluidity. The configuration-constrained PES calcula- 
tions allow us to determine the fission barrier for a given 
configuration. Fig.0|shows the configuration-constrained 
fission barrier for the 7~ isomer in 256 Fm. (It needs to 
be mentioned that an isomer might undergo a more com- 
plex dynamical process during fission, without keeping a 
fixed configuration.) The increase of the fission barrier 
in both height and width implies an increase in the fis- 
sion lifetime. In the calculations of fission barriers, the 
inclusion of the non-axial 7 deformation is important be- 
cause the 7 degree of freedom can significantly affect the 
shape of the fission barrier. In 256 Fm, for example, the 
non-axial PES calculations result in the barrier heights 
being lowered by 2.4 and 2.6 MeV for the g.s. and 7~ 
isomer, respectively, compared with axial calculations. 
Fig. El displays the calculated PES of the g.s. in 256 Fm, 
with the fission trajectory marked. The 7~ isomer has a 
similar fission trajectory. 

Although the production of high-spin states in super- 
heavy nuclei is generally inhibited by their increased fis- 
sion probability, this does not apply when the angular 
momentum is non-collective and aligned with the sym- 
metry axis, as for high-K states. Therefore, there should 
be significant production of high-K isomers in fusion re- 
actions. Two-qp excitations can also appear in odd-mass 
and odd-odd superheavy nuclei, resulting in three- and 
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FIG. 5: Calculated PES for the ground state of 2j6 Fm. The 
dashed line indicates the fission trajectory. The filled circle, 
triangle and square are at the minimum, maximum and saddle 
points, respectively. At each (fh,f) point, the energy has 
been minimized with respect to the /?4 deformation. 

four-qp states, respectively. These have not yet been con- 
sidered in detail. Furthermore, the present investigation 
is limited to deformed superheavy nuclei. However, the 
angular-momentum effects should also apply to spherical 
superheavy nuclei. 

In summary, we have investigated multi-qp states 
in even-even superheavy nuclei using configuration- 
constrained PES calculations. Low-energy high- if de- 



formed configurations are predicted to be present in the 
vicinity of the Fermi surface, resulting in the formation 
of high-spin isomeric states. The high-spin isomerism in- 
creases the barrier to fission and decreases the probabil- 
ity of a-particle emission. For a given class of superheavy 
nuclei with sufficiently short lifetimes, we have discussed 
an inversion of stability, in that the excited states can be 
longer lived than the corresponding ground states. These 
effects result in increased survival probabilities of super- 
heavy nuclei. Therefore, it can be stated that via the 
population of high- if states in experiments, one may be 
able to extend the nuclear chart further into the island of 
superheavy nuclei. The energy calculations of multi-qp 
states and their a decays indicate that the single-particle 
orbits of the deformed Woods-Saxon potential apply also 
to the superheavy mass region. Further studies of high- if 
states will provide essential tests, not only of our specific 
isomer predictions, but also of mean-field models in gen- 
eral. 
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